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I""— ■ Abstract 
On 

f| We present a next-to- leading order calculation for the virtual photoproduction of one 

and two jets in ep collisions. Soft and collinear singularities are extracted using the 



Qh! phase space slicing method. The collinear photon initial state singularity depends 

logarithmically on the mass of the virtual photon and is absorbed into the virtual 



photon structure function. An MS factorization scheme is defined similarly to the 
real photon case. Numerical results are presented for HERA conditions using the 
Snowmass jet definition for inclusive single jet and dijet cross sections. We study 
the dependence of these cross sections on the transverse energies and rapidities of 
the jets. Finally, we compare the ratio of the experimentally defined resolved and 
direct cross sections with recent ZEUS data as a function of the photon virtuality 
P 2 . 
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1 Introduction 



In ep scattering at HERA interactions between photons of small virtuality P 2 and protons 
produce jets of high transverse energy Ex- The presence of this hard scale Ex allows the 
application of perturbative QCD to predict cross sections for the production of two or more 
high-i?T jets which can be confronted with experimental data. This offers the opportunity 
to test QCD and to constrain the structure of the colliding particles. 

At leading order (LO) QCD two distinct processes are responsible for the production 
of jets. In the direct photon process, the photon interacts as a point-like object with 
a parton in the proton, whereas in the resolved process the photon acts as a source of 
partons which then scatter with the partons coming from the proton. Both LO processes 
are characterized by having two outgoing jets of large transverse energy. Studies of dijet 
photoproduction at HERA have shown that both classes of processes are present for the 
case of quasi-real photons, i.e. photons of extremely small virtuality P 2 ~ [I], 0. The 
comparison of theoretical predictions || |j] with these and more recent experimental 
data H have given us some confidence that the parton distribution functions for the real 
photon available in the literature are consistent with the dijet production data. The 
parton content of photons with virtuality P 2 = is reasonably constrained by data from 
deep inelastic scattering J7J. Unfortunately this is not the case for a photon target with 
non-zero, although small, virtuality P 2 . The only measurement for the virtual photon 
structure function available so far has been performed by the PLUTO collaboration at 
PETRA g. They measured the structure function F eff = F 2 + \F L for P 2 < 0.8 GeV 2 
and Q 2 = 5 GeV 2 as a function of x (Q 2 is the virtuality of the probing virtual photon, 
whereas P 2 always denotes the virtuality of the probed virtual target photon). More 
and better data should come from LEP2 f||]. On the theoretical side several models 
exist for describing the Q 2 -evolution equations of the parton distributions and the input 
distributions at Q with changing P 2 [[K| [□], [12| . These constructions use essentially the 



same methods as have been applied for the parton distributions of real photons. This 
way some smooth behavior towards P 2 = 0, where previous results for the real photon 
should hold, is guaranteed. In [|10| this construction allows a calculation of the parton 
distribution functions (PDF) for virtual photons in LO and NLO. It incorporates a purely 
perturbative contribution and a non-pointlike hadronic contribution. Unfortunately these 
PDF's for P 2 7^ are not available in a form that parametrizes the Q 2 evolution. Such 
parametrized PDF's for virtual photons have been presented recently by two groups [O, 



|12| , but unfortunately only in LO. These two models differ somewhat in their method 
of extrapolation to P 2 ^ 0, in the choice of the input distribution and the choice of 
the input scale Q . Furthermore the PDF's of Gliick, Reya and Stratmann (GRS) [0 
present distributions only for Nj = 3 flavors, whereas Schuler and Sjostrand (SaS) give 
PDF's for Nf — 4 flavors [|11|]. We expect these two structure function sets to be detailed 
enough so they can be tested in photoproduction experiments with virtual photons. First 
preliminary data have been presented by the ZEUS |TJ[ and the HI [|TjJ collaborations. 
ZEUS studied the dijet cross section for E T > 4 GeV in the range of < x 7 < 1. Events 
with x 7 > 0.75 are assumed to be dominated by the direct process whereas events with 
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< x 7 < 0.75 give the resolved component. Then the ratio of the direct-enriched and 
resolved-enriched cross section is measured as a function of P 2 for < P 2 < 0.55 GeV 2 . 
In the HI experiment the inclusive one-jet cross section is measured as a function of 
Et and rapidity 77 for various P 2 bins [16|. It is expected that the resolved component 



decreases relative to the contribution from the direct photon processes as the virtuality of 
the photon increases. Some theoretical studies of the single inclusive and dijet inclusive 
cross section in LO have been presented recently [12 , [13 , |l4 . 



It is well-known that in NLO calculations the distinction between direct and resolved 
photoproduction becomes ambiguous. In this order, a large contribution in the NLO di- 
rect cross section is subtracted from the direct component and combined with the LO 
resolved term thus producing the scale (Q = M 7 ) dependence of the PDF's of the pho- 
ton. Therefore both components are related to each other through the factorization scale 
M-y at the photon leg which determines the part of the NLO direct contributions to be 
absorbed into the resolved component. The M 7 dependence of the remaining NLO direct 
contribution cancels to a large extent against the dependence in the resolved cross sec- 
tion coming from the photon structure function. This connection has been worked out 
in detail JlT], |J and studied numerically for real photoproduction with P 2 = fT8| . It 
is clear that this cancellation of the scale dependence must take place also for P 2 7^ 0. 
So, for a consistent calculation up to NLO one needs to superimpose both contributions, 
the NLO direct and at least the LO resolved cross section, both computed with the same 
scale. The subtraction of the large contribution in the direct cross section which is shifted 
to the resolved term is defined only up to finite, non-singular terms in the limit P 2 — > 0. 
These finite terms may be fixed in a way that a smooth behavior towards the limiting 
case of real photons (P 2 = 0) is achieved where these finite terms are usually defined in 
the MS subtraction scheme. For a complete NLO calculation we must include the NLO 
hard scattering parton-parton cross sections for the direct (here one of the ingoing par- 
tons is the photon with virtuality P 2 7^ 0) and for the resolved process together with the 
two-loop evolved structure function of the proton and photon with virtuality P 2 7^ 0. 

In this paper we shall work out the subtraction in the NLO direct contribution and 
superimpose the remaining direct contribution and the resolved cross section up to NLO. 
We study this cross section and the two contributions as a function of P 2 for various 
inclusive one- and two-jet cross sections. Of particular interest is the behavior of the 
resolved component as a function of P 2 and the question at which P 2 the sum of direct 
and resolved cross sections approaches the unsubtracted direct cross section which one 
expects to give the correct description at sufficiently large P 2 . 

The outline is as follows. In section 2 we describe how to subtract the singular terms 
in the NLO direct cross section connected with the collinear singularity of the 7 — > qq 
contribution. Here we also define the finite terms which depend on the subtraction scheme 
used for P 2 = 0. Our result concerning the P 2 dependence of various one- and two- 
jet cross sections are presented in section 3. In this section we also compare with the 
preliminary data from ZEUS Section 4 contains a summary and an outlook for 

future investigations. 
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2 Subtraction of Photon Initial State Singularities 



The NLO corrections to the direct process become singular in the limit P 2 = 0. For 
real photoproduction these photon initial state singularities are usually evaluated, i.e. 
regularized, with the dimensional regularization method in which they are finite for e = 
(4 — d)/2. The singular contributions appear as poles in e multiplied by the splitting 
function P 9i< - 7 |19j . These initial state singularities are absorbed into the PDF of the real 
photon (P 2 = 0). With no further finite terms subtracted (for e — > 0) this defines the MS 
factorization scheme which must be consistently applied also for the NLO evolution of the 
photon PDF. For P 2 ^ the corresponding contributions appear as terms proportional 
to ln(P 2 /s), yfs being the cm. energy of the photon-parton subprocess. These terms are 
finite as long as P 2 ^ and can be calculated with d — 4 dimensions. Since for small 
P 2 these terms are large they are absorbed as in the case P 2 = into the PDF of the 
virtual photon which is present in the resolved cross section. Concerning finite terms (for 
P 2 —>■ 0) which may also be subtracted together with the singular terms we have the same 
freedom as in the case P 2 = 0. We fix these finite terms such that they agree with the MS 
factorization in the real photon case. To achieve this we must compare the contributions 
originating from the photon initial state singularities in the two cases P 2 = and P 2 ^ 0, 
which we shall do in the following. 

For this purpose we must isolate the singular terms from the photon initial state 
when the photon is collinear with one of the outgoing quarks or antiquarks. The relevant 
processes are labeled K 1 ,K 2 ,K 3 and as specified in Tab. 1. We do not separate the 
contributions according to color factors as in |19|| . To make the comparison with [[0| 



possible we have given the contribution to the relevant processes also in the notation used 
|i~9| . Using the same definitions of momenta and variables as in [[Tj| the integration of 



m 



the 2 — > 3 matrix elements over the phase space dPS< p > yields for P 2 = the following 
result (i = 1,2,3,4): 



dZn 



J dPS {r] H Ki = J ^eV/i 4e 8(l - e 




Table 1: Classification of 2 — ► 3 matrix elements 



Class 


Process 


Class in [19] 


Ki 


iq -> qgg 


h + h 


K 2 


19 -> qqg 


h + h 


K 3 


79 -> qqq 


h + h + h 


K 4 


7g — > qq'q' 


h 
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The result for the K, L can easily be read from the results in appendix C of |T9j by adding 
the corresponding sums of ij's according to Tab. 1. The result is written in the following 
form 



K 2 



K A 



= 2C F MU 1 (s,t,u) 

= -2M Ui (t, s, u) + U x (u, s, t) 

= iC F M U a (s,t,u) + U 2 (t,s,u) 

+ ±(1/3(3, t, u) + cycl. permutations in s, t, u) 



2C F M 



U 3 (s, t, u) + cycl. permutations in s, t, u 



where M is defined as 



M 



1 1 



2N, 



-Pn. 



c 



2N, 



-Pn 



c 



- 7 (z a ) In 



[1 - Zg 



+ f 



In (6) z a = e [x a , 1] and the splitting function 



(2) 
(3) 

(4) 
(5) 

(6) 



P qi ^(z a ) = 2NcQi 



(7) 



The functions Ui(s,t,u) are the LO parton-parton scattering cross sections related to 
the various processes as shown in Tab. 2. Processes, which are related by crossing are 



omitted, the complete list is given in [IS], q and q' denote quarks with different flavors. 
The explicit expressions for the Ui and their dependence on color factors are given in 
||19|| . The factor M contains the characteristic singularity proportional to 1/e which is 
subtracted by absorbing 



R, 



Z a ,M* 



-Pn. 



r(i-e) 
r(i - 2 £ ) 



'47T/i 2 



into the PDF of the photon F a / y (x a , M^) (see |I9| ). This subtraction produces the fac- 
torization scale (M 7 ) dependence of the photon PDF and gives the finite contributions to 
the cross section which are given by (2)-(5) with M replaced by Mjf§: 



AI- 



MS 



2N< 



c 



>a) In 



K y s s(i 



+ 



01 



(9) 



The variable y s which appears in (6) and (9) is the invariant mass cut-off used to isolate 
the collinear singularity in the 7 — > qq splitting. 

The procedure for virtual photons with virtuality P 2 is completely analogous. First 
one calculates the 2 — > 3 matrix elements, but now with P 2 7^ and decomposes them into 
terms with the characteristic denominator from 7 — > qq splitting which become singular 
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Table 2: LO parton-parton scattering matrix elements 



Process 


Matrix elements |.M| 2 = SNcCpg^Ui 


QQ ~^ 99 


C/i(s,t,it) 


qq' -> qq' 


U 2 (s,t,u) 


QQ -> 99 


U 2 (s, t, u) + C/ 2 (s, it, £) + C/ 3 (s, i, w) 



in the limit P 2 — > 0. These terms after phase space integration up to a cut-off y s as in (1) 
have the same structure as (2)-(5) with the LO parton-parton matrix elements factored 
out. The integration can be done with e = since P 2 ^ 0. The phase space in (1) then 
contains the additional P 2 -dependent factor / 

P 2 (l - z ) 

which reduces to 1 for P 2 = 0. The factor M in the equations (2)-(5) takes the simple 
form 

M=^(l + ^)P^a) (11) 



2N C V z a P 2 , 

which is singular for P 2 = as to be expected. This singularity is absorbed into the PDF 
of the virtual photon with virtuality P 2 . Instead of (8) the subtraction term is: 



/ M 2 \ 

M ') = ln [ p2 {1 l Za) ) p «-rW - N cQi ■ (12) 
After this subtraction the remaining finite term (for P 2 — > 0) in M yields 

M^-^M^ ^^.J .f . d3) 

In addition to the singular term \n(M 2 /P 2 ) we have subtracted in (12) two finite terms in 
order to achieve in (13) the same result as in (9) for P 2 = 0. Therefore we call this form 
of factorization the MS factorization for P 2 ^ 0. It is defined by the requirement that the 
remaining finite term is equal to the finite term in (2)-(5) with M replaced by Mjfg as 
was obtained in the calculation for real photons. With this definition of factorization in 
the case P 2 ^ we make sure that we obtain the same NLO corrections as in where 
P 2 = 0, when in the P 2 ^ calculation we choose P 2 extremely small. So, for the actual 
calculations we apply the formulas (2)-(5) where M is now given by (13). We note that 
for averaging over the spin of initial gluons and photons we apply a factor 1/(2 — 2e), that 
gives rise to some additional finite terms which should be included in (9) when expanded 
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in e. We have taken these terms into account also in our calculation. This completes the 
calculation of the contribution from the photon initial state singularity. 

It is clear that in NLO the PDF for the virtual photon must be given also in the MS 
factorization scheme. In ref. |H| the PDF is constructed in the so-called DIS 7 scheme, 



which is defined as for real photons (P 2 = 
MS scheme PDF in the following way fLO 



0). This distribution function is related to the 



F a 



th (x, M 2 ) DISl = F ah (x, M 2 )^ + 6F 



^ /7 (x,M 7 2 ) 



(14) 



where for a = q^q\ or g: 
5F qih {x,M 2 ) = 

6F gh (x,M 2 ) - 



-N c 



1 



6Fq. h (x,M, 
a 





P ? ^ 7 (x)ln 



Q 2 8x{l-x)-Q 2 



(15) 



If the PDF in this scheme is used to calculate the resolved cross section one must transform 
the NLO finite terms in the direct cross section. This produces a shift of M(P 2 )jjg 
given in (13) by the same expression as (15). Then the relation is: 



as 



M, 



!>IS~,(P Z 



M(P 2 )^ - N< 



N c Pq ^ Za 



In 



1 - z a 



Q?8* (l 



Ql 



• (16) 



All other singular terms in the real corrections, i.e. final state singularities and the 
contributions from parton initial state singularities have been calculated by Graudenz in 
connection with the NLO corrections for jet production in deep inelastic ep scattering. 



They can be taken from his work |2(J together with the virtual corrections up to 0{aa s ). 
(For related work see pip .) All these contributions are calculated in dimensional regular- 
ization. The appearing singularities in 1/e cancel when all singular terms are added. The 
remaining finite terms enter the NLO corrections for the jet cross sections. These finite 
terms depend on the phase space slicing parameter y s which is introduced to separate the 
singular regions of final and initial state infrared and collinear divergences. 



3 Inclusive One- and Two-Jet Cross Sections 

In this section, we present some characteristic numerical results for one- and two-jet cross 
sections as a function of the virtuality P 2 . We consider the contributions of the direct and 
resolved components and their sum. Partly we shall compare the NLO results with LO 
predictions. The input for our calculation is the following. We have chosen the CTEQ3M 



proton structure function |22| which is a NLO parametrization with MS factorization and 



Atth = 239 MeV. This A value is also used to calculate a s from the two-loop formula at 
the scale fi = Et- The factorization scales are also M 7 = M p = Et- We also need the 
parton distribution of the virtual photon F a /^. For this we choose either the GRS [|l^] set 



or the SaSIM set of Schuler and Sjostrand ||11|| . Both sets are given in parametrized form 
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for all scales M 2 so that they can be applied without repeating the computation of the 
evolution. Unfortunately both sets are given only in LO, i.e. the boundary conditions for 
P 2 = M 2 and the evolution equations are in LO. In [ 10 1 PDF's for virtual photons have 
been constructed in LO and NLO. Distinct differences occur for larger P 2 and x > 1CT 3 
which is mainly due to the different NLO perturbative boundary condition at P 2 = M 2 , 
which does not exist for the real (P 2 = 0) photon structure function. Since neither 
of the two PDF's is constrained by empirical data from scattering on a virtual photon 
target we consider these LO distribution functions as sufficient for our exploratory studies 
on jet production and treat them as if they were obtained in NLO. In particular, we 
shall use the transformation formulae (14) for going from the DIS 7 to the MS-scheme 
PDF, which makes sense only in NLO. Then the parametrization 0] is considered as the 
parametrization in the DIS 7 scheme. In the PDF of GRS the input scale is Qo = 0.5 GeV 
and the restriction P 2 < Q 2 /5 is implemented as to fulfill the condition A 2 < ? 2 < Q 2 . 
In addition the PDF of GRS is constructed only for Nf = 3 flavors. The production of 
the heavier c and b quarks is supposed to be added as predicted by perturbation theory 
of fixed order with no active c and b quarks in the proton and photon PDF's. In LO this 
amounts to adding the processes i*g ~^ °c an d 1*9 —> bb to the cross section, keeping 
m c , nib ^ 0. Since in this work we are primarily interested in studying the sum of the 
direct and resolved contributions and the influence of the consistent subtractions of the 
NLO direct part we refrain from adding the LO or NLO cross sections for direct heavy 
quark production as suggested in JR], |T2|. So, our investigations in connection with the 
GRS parametrization of the virtual photon PDF are for a model with three flavors only. 
For consistency we take also Nf = 3 in the NLO corrections and in the two-loop formula 
for a s . To overcome this problem we studied the relevant cross sections also with the 
virtual photon PDF's of [TT|. Here the c quark is included as a massless flavor in the PDF 
which undergoes the usual evolution as the other massless quarks except for a shift of the 
starting scale Qo- This Nf = 4 PDF is considered only in the parametrization SaSIM 
with Qo = 0.6 GeV that is in the MS scheme. 

The cross sections we have computed are essentially for kinematical conditions as 
in the HERA experiments. There, positrons of E e = 27.5 GeV produce photons with 
virtuality P 2 which then collide with a proton of E p = 820 GeV. The momentum transfer 
to the photon is q = p e — p' e with P 2 = (— q 2 ). The energy spectrum of the virtual photons 
is approximated by 

dEygy) _ a l + (l-yf 1 
dP 2 2tt y P 2 
with 

V=**-*^ (18) 

PPe E e 

being the fraction of the electron energy transferred to the photon, when the virtuality 
P 2 <C E 2 , . The momentum of the incoming proton is p. The approximation for the virtual 
photon spectrum is used for the calculation of the resolved and the direct cross section. 
The expression ( |j~7|) factorizes in the cross section for ep — ► e'X with arbitrary final state 
X if one neglects the longitudinal virtual photon terms. After integration over P 2 between 



(17) 
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Pmin < P 2 < Pmax with P^ in := m 2 y 2 /{ 1 — y) , where m e is the electron mass, one obtains 
the Weizsacker- Williams formula as used for calculations with untagged electrons as in 
Jl|. [TjJ, where P 2 ~ dominates. The cross section for the process ep —>■ e'X is then given 
by the convolution 

dajep -> e'X) T, dF l/e {y) , 

Vinin 

where da^p — > X) denotes the cross section for 7*p — > X with transversely polarized 
photons of energy go — E e y — P 2 / (4E e y) , if the transverse component of q is neglected. 
To have the equivalent conditions as in the ZEUS analysis we choose y min = 0.2 and 
Umax = 0.8. In the computation of the resolved cross section the approximation q = E e y 
is inserted for the energy of the virtual photon, whereas for the direct photon cross section 
the exact relation for go is taken into account through the kinematic relations. 

All further calculations proceed in the following way. For both, the direct and the 
resolved cross section, we have a set of two-body contributions and a set of three-body 
contributions. Each set is completely finite, as all singularities have been canceled or 
absorbed into PDF's. Each part depends separately on the phase space slicing parameter 
y s . The analytic calculations are valid only for very small y s , since terms 0(y s ) have been 
neglected in the analytic integrations. For very small y s , the two separate pieces have 
no physical meaning. In this case the (ln7/ s )-terms force the two-body contributions to 
become negative, whereas the three-body cross sections are large and positive. In || we 
have plotted such a cross sections for direct real photoproduction at y s = 10~ 3 . When 
the two-body and three-body contributions are superimposed to yield a suitable inclusive 
cross section, as for example the inclusive one- or two-jet cross section, the dependence 
on the cut-off y s will cancel. The separation of the two contributions with the slicing 
parameter y s is a purely technical device in order to distinguish the phase space regions, 
where the integrations are done analytically, from those, where they are done numerically. 
We have checked, by varying y s between 10~ 4 and 10~ 3 , that the superimposed two- and 
three-body contributions are independent of y s for the inclusive single- and dijet cross 
sections. 

First, we study the inclusive single-jet cross section. To calculate it we must choose a 
jet definition, which combines two nearly collinear partons. We adopt the jet definition 



of the Snowmass meeting p3| . According to this definition, two partons i and j are 



recombined, if Rij < R, where Rij = y (rji — r/j) 2 + (0j — 0j) 2 and 77 j, <pj are the rapidity 
and the azimuthal angle of the combined jet respectively, defined as 

Etj = Et x + Et 2 (20) 

nj = E ^ + E ^ 2 , (2i) 



Etj 
E Tl 4>i + E T2 ( 



Ej- 



(22) 



j 



and R is chosen as in the experimental analysis. So, two partons are considered as two 
separate jets or as a single jet depending whether they lie outside or inside the cone with 
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radius R around the jet momentum. In NLO, the final state may consist of two or three 
jets. The three-jet sample contains all three-body contributions, which do not fulfill the 
cone condition. The cone constraint is evaluated in the HERA laboratory system as for 
real photoproduction (P 2 = 0) and in the experimental analysis. The calculation of the 
resolved cross section proceeds as for real photoproduction, i.e. the transverse momentum 
(qr) of the virtual photon is neglected so that the virtual photon momentum is in the 
direction of the incoming electron and go = E e y. The cross section for direct photons, 
in which the initial state singularity is subtracted, as specified in the previous section, is 
given in the center-of-mass system p + q = and transformed into the HERA laboratory 
system, again neglecting and other small terms proportional to P 2 . Then the relation 
between the rapidity r] cm of the jet in the c.m.-system and laboratory system is as for real 
photoproduction, which is 

1 E 

V = Vcm + 77 In ■ (23) 

2 yE e 

This 7] and the corresponding azimuthal angle of the partons in the final state is also used 
for evaluating the jet definition given above. 

In Fig. la, b, c, the results for d 2, a / dETdr/dP 2 are shown integrated over r\ in the 
interval —1.125 < f] < 1.875 (these boundaries are employed in the ZEUS analysis fL5f ) 
as a function of Ex for the three values of P 2 = 0.058, 0.5 and 1 GeV 2 . The smallest value 
of P 2 has been chosen in such a way that it reproduces the P 2 ~ results. Inserting 
P 2 = 0.058 GeV 2 into the unintegrated Weizsacker- Williams formula corresponds to the 
one integrated in the region P^ in < P 2 < P^ax = 4 GeV 2 . For all three P 2 the cross 
section is dominated by the resolved component at small Et- Near Et = 20 GeV the 
direct component, which is the direct cross section without the subtraction term, denoted 
by Dir s , is of the same magnitude as the resolved cross section. As a function of P 2 the 
cross section d 3 a / 'dE^d^dP 2 decreases more or less uniformly in the considered Et range 
with increasing P 2 . One can see that this decrease is stronger for small Et as compared to 
large Et- In these and the following results the cone radius is R = 1. The corresponding 
rapidity distribution for a fixed Et will be shown later in Fig. 5a, b, c, d for the SaSIM 
photon PDF. 

In Fig. 2a, b, c we studied the r] distribution of the Dir s contribution at fixed Et = 7 
GeV and the same P 2 values as in Fig. 1 in comparison with two approximations, namely 
the LO cross section and the NLO cross section with P 2 ^ only in the photon flux 
equation (17) and the rest of the cross section evaluated at P 2 = as in [|nj. As to 
be expected this approximation is very good for P 2 = 0.058 GeV 2 . At the larger P 2 
however it overestimates the cross section and should not be used. This means that the 
P 2 dependence of the direct part, although the strongest logarithmic P 2 dependence has 
been subtracted, should be taken into account. Of course, in the sum of the resolved 
and the direct cross section the difference is small as long as the resolved part dominates 
which is true for the smaller EtS. The LO prediction, which is evaluated with the same 
structure functions and a s value as the NLO result, only the hard scattering cross section is 
calculated in LO, is smaller than the NLO result, the difference decreasing with increasing 
P 2 . Of course, this finding depends on the chosen value of R. The NLO cross section 



11 



depends on R, whereas the LO curve does not. Therefore estimates of the inclusive cross 
section with LO calculations can be trusted only for large cone radii. 

It is clear that the resolved and the direct cross sections decrease with increasing P 2 
for fixed rj and Et- It is of interest to know how the ratio of resolved and the direct cross 
section behaves as a function of P 2 . This is a well defined quantity in LO. The variation 
of this ratio with P 2 up to P 2 = 4 GeV 2 for E T = 7 GeV and r\ = 2, 1, and -1 is plotted 
in Fig. 3a, b, c, d. As we expected this ratio decreases most strongly for t] = 2, since 
in the r] > region the resolved component dominates whereas the direct cross section 
peaks for negative 77's (see Fig. 2a, b, c). If we decrease r] towards the backward direction 
this ratio diminishes more or less uniformly for all P 2 . In NLO this ratio depends on the 
scheme chosen for the photon PDF. In the DIS 7 subtraction scheme terms in the photon 
PDF are shifted to the direct cross section as follows from (14)-(16) in section 2. This 
necessarily changes the ratio of the resolved to the direct cross section, not only in the 
absolute value but also in the dependence on P 2 . For all r/'s this ratio is quite different 
from its LO value. The difference between the MS and DIS 7 scheme is small for 77 = 0, 1 
moderate for rj = 2 and of the order of 1.5 at rj = — 1. Except at rj = —1 the ratio is 
always larger in the DIS 7 scheme than in the MS scheme. Apart from the fact that the 
ratios plotted in Fig. 3a, b, c, d cannot be measured directly, they are scheme dependent 
in NLO and have very large corrections when going from LO to NLO. In Fig. 3a, b, c 
we have also plotted the ratio for the case that the P 2 dependence is neglected in the 
cross section 7*p — > X and taken into account only in the photon propagator (denoted 
NLO(P 2 = 0)). This approximation gives a result very similar to the LO curve showing 
that the NLO corrections are more important for P 2 7^ 0. Of course at P 2 — > this 
approximation is equal to the NLO result in the MS scheme. 

The results shown so far are for a model with three flavors only and therefore should 
not be compared to the experimental data except when the contribution from the charm 
quark is added at least in LO. A more realistic approach is to use the photon PDF's 
SaSIM of ref. [|TIJ which are constructed for four flavors. The results for d 3 a/dETdi]dP 2 
integrated over 1] e [-1.125, 1.875] as a function of E T for P 2 = 0.058,0.5 and 1.0 GeV 2 
are presented in Fig. 4a, b, c. These curves can be compared with the results in Fig. la, 
b, c obtained with the GRS photon distribution with Nf = 3. The results for the sum of 
the resolved and direct contributions change between 10% and 30% in the small Ej- region 
and approximately 50% in the large Et region. As one can see the larger cross section for 
Nf = 4 results primarily from the direct component. Since the direct component is more 
important for large Et than for small Et the increase is stronger in the former region. 

Of interest are also the rapidity distributions for fixed Et- These are shown for Et = 7 
GeV as a function of 77 between — 1 < 77 < 2 choosing P 2 = 0.058, 1,5 and 9 GeV 2 in 
Fig. 5a, b, c, d. We show the subtracted direct cross section and the resolved cross 
section and their sum for the photon PDF as in Fig. 4a, b, c. Expectedly the resolved 
component has its maximum shifted to positive 77's as compared to the direct component. 
The direct component falls off with increasing 77 quite strongly. This comes from the 
subtraction of the (InP 2 ) terms as is apparent when we compare with the unsubtracted 
direct cross section, labeled "Dir" in Fig. 5a, b, c, d. The sum of the resolved and the 
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direct (subtracted) cross section is more or less constant for the smaller P 2 values and 
decreases with increasing 77 for P 2 = 5 and 9 GeV 2 . 

For large enough P 2 we expect the unsubtracted direct cross section to be the correct 
one. In this region the subtraction term (12) must approximate the PDF of the photon 
rather well. We have checked this by a direct comparison. With increasing P 2 the full 
direct cross section (Dir) approaches the Sum=Res+Dir s better and better. But even at 
P 2 = 9 GeV 2 the two cross sections differ at 77 = 2 still by approximately 30 %. However, 
in the backward direction at 77 = —1 we also see a difference. In this region, which 
corresponds in the photon PDF to the region x 1 ~ 1, where the perturbative component 
dominates, we do not expect a deviation. This may be caused by our approximation of 
neglecting the transverse momentum of the virtual photon in the relation (22) and in 
the calculation of the resolved cross section. 

In order to compare with preliminary data for the dijet cross section presented by the 
ZEUS collaboration |L5] we calculated the inclusive dijet cross section d 4 a / dETdr/id^dP 2 
as a function of P 2 . Here Et is the transverse energy of the measured or trigger jet with 
rapidity 771. 772 denotes the rapidity of another jet such that in the three-jet sample these 
two jets have the largest Et, i.e. Et 17 Et 2 > Et 3 - The calculational procedure is analogous 
to the inclusive single-jet cross section and follows closely the work presented in [|J for 
P 2 = 0. Since inclusive two-jet cross sections depend on one more variable they constitute 
a much more stringent test of QCD predictions than inclusive one-jet cross sections. We 
could predict distributions in rji and 772 for fixed Et or distributions in Et for various 
values of the two rapidities rji and 772 in the same way as in [|J for P 2 = 0. Since no such 
information is expected from the experiment in the near future we calculated only the Et 
distribution with the two rapidities integrated over the interval —1.125 < 771,772 < 1.875 
following the constraints in the ZEUS analysis [15]. The results for P 2 = 0.058,0.5 and 
1.0 GeV 2 are shown in Fig. 6a, b, c, where the full curve is d A a j dETdrjid^dP 2 as a 
function of Et integrated over 7/1 and 7/2 in the specified interval and for 0.2 < y < 0.8. 
The functional dependence on Et does not change as a function of P 2 , only the absolute 
value of the cross section decreases with increasing P 2 . 

In Fig. 6a, b, c we show also the cross sections for enriched direct and resolved 7 
samples. These cross sections are labeled "Dir" (dashed curve) and "Res" (dotted curve), 
respectively. The two cross sections are defined with a cut on the variable x° bs where x obs 
is defined by 

x° bs = ^ T ; . (24) 



This variable measures the fraction of the proton energy that goes into the production of 
the two hardest jets. The "Dir" curve gives the cross section for x° bs > 0.75. This cut on 
x° bs leads to an enrichment of the direct component of the cross section, since exclusive 
two-jet events from the direct process have x° bs = 1. The curve labeled "Res" gives the 
cross section in the complementary region x° < 0.75 which characterizes the enriched 
resolved 7 sample. The sum of the Dir and Res curves is equal to the full cross section 
d A o 1 1 dETdr\\dr\2dP 2 with no cut on x obs . It must be emphasized that both cross sections, 
whether x obs > 0.75 or x° bs < 0.75, contain contributions from the direct and resolved 
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part. Actually the resolved contribution in the enriched direct sample (x° bs > 0.75) 
is non-negligible. The results in Fig. 6a, b, c are for the GRS parton distributions in 
the MS scheme. As to be expected, with increasing P 2 the full cross section is more 
and more dominated by the Dir component, in particular at the larger Et- This means 
that the cross section in x° bs < 0.75 decreases stronger with P 2 than in the x° bs > 0.75 
region. This could be studied experimentally by measuring the ratio of the two cross 
sections as a function of P 2 for fixed Et- This has not been done yet. Instead, the 
ZEUS collaboration [15] measured the ratio r =Res/Dir, where Res and Dir are the cross 
sections as defined above, but integrated over E Tl ,E T2 > 4 GeV. With this integration 
cut on the transverse momenta of the two hardest jets, the transverse momentum of the 
third jet may vanish which affords particular constraints on the remnant jets. This is 
treated as in [Q]. Furthermore we replaced the GRS photon PDF by the SaSIM photon 
PDF which is for Nf = 4 flavors. With these assumptions we calculated the ratio r as a 
function of P 2 up to P 2 = 0.6 GeV 2 and compared it with the ZEUS [|T5J data in Fig. 7 
in LO (dotted curve) and NLO (full curve). The theoretical NLO curve agrees quite well 
with the data at P 2 > 0.25 GeV 2 but not with the measured point at P 2 ~ 0.2 GeV 2 
and at P 2 = 0.058 GeV 2 corresponding to the photoproduction case. Of course the ratio 
r for P 2 ~ is much more precise and lies 30 % higher than the predicted cross section. 
This disagreement at P 2 ~ is to be expected since at this value of P 2 the measured 
inclusive dijet cross section for the enriched resolved 7 sample is larger than the predicted 
cross section for a small Et cut M. As in M we attribute this difference between theory 
and experimental data to additional contributions due to multiple interactions with the 
proton remnant jet in the resolved cross section not accounted for by our NLO predictions. 
This underlying event contribution is reduced with increasing E T im and/or smaller cone 
radii R < 1. As it seems, for larger P 2 , the underlying event contribution is also reduced. 
Of course this could be studied more directly by measuring rapidity distributions for the 
enriched resolved 7 sample as was done at P 2 ~ in [Q. 



4 Conclusions 

We have calculated inclusive single jet and dijet cross sections for photoproduction with 
virtual photons. The direct and resolved contributions were calculated in next-to-leading 
order QCD using the phase-space slicing method. They were folded with the uninte- 
grated Weizsacker- Williams approximation and existing LO parametrizations for the vir- 
tual photon parton densities. The collinear singularity in the direct photon initial state 
was integrated analytically up to an invariant mass cut-off y s . Contrary to real photons, 
this specific singularity is not regulated in the dimensional regularization scheme but by 
the mass of the virtual photon P 2 leading to a logarithmic dependence on P 2 . This log- 
arithmic term is absorbed into the virtual photon structure function rendering the latter 
scheme and scale dependent. The remaining finite contribution is constructed in such 
a way that the corresponding real photon term is obtained in the limit P 2 = in the 
MS scheme. Similarly to the construction of virtual photon structure functions by GRS 
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and SaS, our calculation of the hard scattering cross section then provides a consistent 
extension from P 2 = to small, but non-zero P 2 . 

We presented y s -cut independent results using the Snowmass jet definition and HERA 
conditions for distributions in the transverse energy and the rapidity of the observed jet 
and in the photon virtuality. For the case of very small P 2 , we found good numerical agree- 
ment with the predictions for real photoproduction. At P 2 = 9 GeV 2 , the unsubtracted 
direct contribution corresponding to the case of deep inelastic scattering approximates the 
sum of the subtracted direct and resolved contributions quite well. A small discrepancy 
remains in the forward region, where the resolved contribution is dominant. As in the case 
of P 2 = 0, the resolved component dominates at low E T and in the forward region. The 
NLO effects were demonstrated to be important in the ratio of resolved and direct cross 
sections as a function of P 2 . Since the theoretical separation between direct and resolved 
is artificial, some scheme dependence remains here unless both contributions are defined 
by an experimental cut on x° bs in dijet cross sections. Then the corresponding ratio shows 
significant NLO effects and good agreement with the available ZEUS data for P 2 > 0.2 
GeV 2 . The disagreement below this value can be attributed to additional contributions 
coming from, e.g., multiple scattering between the photon and proton remnants. 

Future investigations on virtual photoproduction will require more data on single 
inclusive jet production as they exist for P 2 = and at larger transverse energies. With 
luminosity permitting, a detailed dijet analysis of an infrared safe cross section such as 
d 4 a I ' dExdriid^dP 2 , where the transverse energies of the two jets are not cut at exactly 
the same value, will provide a much improved insight into the structure of the virtual 
photon. Furthermore, choosing a fc^-cluster-like jet definition with smaller cone radii 
will reduce both the uncertainties in the jet algorithm and in the underlying event. On 
the theoretical side, one possible improvement is the correct treatment of the transverse 
momentum of the incoming photon for larger P 2 including a correct transformation from 
the photonic cm. frame to the HERA laboratory system. For a consistent NLO treatment, 
the inclusion of NLO parton densities for the photon is necessary. These are, however, 
needed in a parametrized form and should also be studied in correlation with deep inelastic 
cy* scattering data. 
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Figure la: Single-jet inclusive cross section integrated over rj G [—1.125,1.875] for the 
virtuality P 2 = 0.058 GeV 2 . The MS-GRS parametrization with N f = 3 is chosen. The 
solid line gives the sum of the subtracted direct and the resolved term. 
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Figure lb: Same as figure la with P 2 = 0.5 GeV 2 . 
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Figure 2a: Single-jet inclusive cross sections for E? = 7 GeV and P 2 = 0.058 GeV 2 . The 
MS-GRS parametrization with Nf = 3 is chosen. Only the direct part with subtraction 
(Dir s ) is plotted. The solid line gives the LO contribution. The dashed curve is the full 
NLO cross section, whereas the dotted curve gives the NLO cross section, where the NLO 
matrix elements have no P 2 -dependence. 




Figure 2b: Same as figure 2a with P 2 = 0.5 GeV 2 . 
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Figure 3a: The ratio of the resolved to the subtracted direct contribution in LO and NLO 
for the GRS parametrization, both in the MS- and the DIS 7 -scheme for r] = 2 and Et = 7 
GeV. The dashed curve is for the NLO-matrix elements with P 2 = for comparison. 
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Figure 3b: Same as figure 3a with rj 
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Figure 3c: Same as figure 3a with r\ = 0. 
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Figure 3d: Same as figure 3a with rj — — 1. 
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Figure 4a: Single-jet inclusive cross section integrated over rj G [—1.125,1.875] for the 
virtuality P 2 = 0.058 GeV 2 . The MS-SaSIM parametrization with N f = 4 is chosen. The 
solid line gives the sum of the subtracted direct and the resolved term. 
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Figure 4b: Same as figure 4a with P 2 = 0.5 GeV 2 . 
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Figure 5a: Comparisons of single-jet inclusive cross sections for = 7 GeV and the 
virtuality P 2 = 0.058 GeV 2 . The MS-SaSIM parametrization with N f = 4 is chosen. The 
solid line gives the sum of the subtracted direct and the resolved term. The dash dotted 
curve is the direct contribution without subtraction. 




Figure 5b: Same as figure 5a with P 2 = 1 GeV 2 
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Figure 6a: Dijet inclusive cross section integrated over 771,772 £ [—1.125,1.875] for the 
virtuality P 2 = 0.058 GeV 2 . The MS-GRS parametrization with Nf = 3 is chosen. The 
solid line is the sum of the direct and the resolved contribution. The dashed line is the 
direct-enriched contribution with x° bs > 0.75 and the dotted curve is the resolved enriched 
contribution with x° bs < 0.75. 
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Figure 6b: Same as figure 6a with P 2 = 0.5 GeV 2 . 
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Figure 6c: Same as figure 6a with P 2 = 1.0 GeV 2 . 
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Figure 7: The ratio of the resolved-enriched to the direct-enriched contributions as cal- 
culated in Fig. 6a, b, c, integrated over E Tl ,E T2 > 4 GeV in LO (dotted) and NLO (full) 
for the SaSIM parametrization with Nf — 4 compared with ZEUS data. 
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